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Abstract 

If one randomly punches holes in an infinite tensed membrane, when does the tension 
cease to exist? This problem was introduced by R. Connelly in connection with applications 
of rigidity theory to natural sciences. We outline a mathematical theory of tension based 
on graph rigidity theory and introduce several probabilistic models for this problem. We 
show that if the "centers" of the holes are distributed in according to Poisson law with 
parameter A > 0, and the distribution of sizes of the holes is independent of the distribution 
of their centers, the tension vanishes on all of for any value of A. In fact, it follows from 
a more general result on the behavior of iterative convex hulls of connected subsets of 
when the initial configuration of subsets is distributed according to Poisson law and the sizes 
of the elements of the original configuration are independent of this Poisson distribution. 
For the latter problem we establish the existence of a critical threshold in terms of the 
number of iterative convex hull operations required for covering all of W^. The processes 
described in the paper are somewhat related to bootstrap and rigidity percolation models. 
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1 Introduction 



Let M be a tensed membrane (film) clamped on the boundary. If one makes a small convex hole 
in this membrane, the tension redistributes over the rest of the membrane. However the situation 
is more complicated when there are several hole in the membrane. It is rather clear that if you 
have a non-convex hole or a couple of convex overlapping holes, tension ought to vanish on the 
convex hull of this set (see Fig. 1.1). 




Figure 1.1 

However, it is less intuitive that tension may vanish at some subset of the complement of a 
collection of convex non-overlapping holes. 




Figure 1.2 



For example, the convex hull of three holes shown on Fig. 1.2 cannot support tension which 
may be verified in practice with a sheet of some elastic material and scissors. Therefore, if 
the area where tension ought to vanish is interpreted as defective, all three polygons on Fig. 
1.1 ought to coalesce into one big defect. Such coalescence effect of a "pinwheel configuration" 
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was first noticed by Connelly (1998). A mathematical explanation of this fact will be given in 
Subsection 1.2. 

In the following subsection we give a formal treatment of tension in the language of graph 
rigidity. Our mathematical treatment of a tensed medium is discrete in nature and based on the 
ideas from rigidity theory, as opposed to classical methods of mechanics of continuous media that 
use smooth functions, tensor fields, and differential equations for describing elastic properties of 
materials. Notice that in the planar case holes with smooth boundary can be well approximated 
by polygons and that the theory based on graph rigidity seem is better suited for computer simu- 
lation. The existence of a "discrete" tension in the complement of the holes has some important 
implications for rigidity properties of large planar graphs realized in the complement of the holes. 
This graph theoretic treatment of tension was introduced in (Connelly, Mitchel, Rybnikov (1998, 
1999)). In this paper we are primarily interested in the behavior of the system with very large 
or infinite number of holes. More specifically, we look at the question of existence of a subset 
supporting tension in the complement of an infinite set of holes in . It follows directly from 
definitions (see Subsection 1.1) that if such set exists, it is infinite and a.s. unique, which is, 
apparently in good agreement with intuition. We consider a few models for the distribution of 
holes in an infinite membrane (Section 2), and completely resolve the question of the existence 
of a tensed subset for a "naive" Poisson model suggested by Connelly. In this model holes are 
distributed according to a Poisson law, and their shapes are IID. In fact we prove a somewhat 
stronger statement, which we hope may have some reasonable interpretation (see Subsection 1.1). 
The percolation models that we describe in this paper are closely related to so-called "bootstrap 
percolation" introduced on trees by Chalupa, Leath and Reich (1979) and later on d-dimensional 
lattices by Kogut and Leath (1981). In these models, points are independently occupied with a 
low density and the resulting configuration is taken as the initial state for dynamics based on 
some collection of local rules, in which the occupation status of a point is updated according to 
the configuration of its neighbors. The rigorous analysis of this model have been conducted by 
Aizenman and Lebowitz (1988). Another work, showing that the critical probability Pc = 1 for 
some bootstrap percolation models was by van Enter (1987). Since then, there has been many 
publications devoted to this topic. For the latest results on bootstrap percolation see Dehghan- 
pour and Schonmann (1997) and Gravner and McDonald (1997). Gravner and Schonmann (1999, 
personal communication) had some ideas how to extend their results for continuous case, but as 
far as the authors know, there has been no published results on the model which we describe 
and analyze below. 

1.1 Rigidity and its Applications 

Over the past two decades 2-dimensional and 3-dimensional random central-force networks have 
been used by physicists for modeling the elastic behavior of glasses within the framework of effec- 
tive medium theory (see Duxbury, Thorpe. Jacobs, Moukarzel (1983, 1995, 1996, 1998, 1999)). 
It turns out that real glasses are well represented by generic random networks. The success of 
these methods resulted in good characterization of elastic properties of glasses like Ge^AsySei^y 
(Thorpe (1983)). The rigidity analysis of central-force networks has been also used for charac- 
terization of physical properties of other substances such has proteins and semiconductors (see 
Thorpe and Duxbury (1999)). 

While discussing the property of a membrane to be able to support tension, we keep in mind 
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(large) networks in the complement of the holes. It turns out that the existence of tension in 
the compliment of the holes implies some important properties for a network realized in the 
complement. At this point we need to introduce some mathematical terminology. Notice that in 
mathematics of rigidity there is a tendency to use term framework instead of network preferred 
by physicists. 

A bar-and-joint framework is a realization of a graph in W'' . Denote by F{E; V, Vq) a frame- 
work in R^'with the edge set E, and the vertex set V, with pinned (fixed in W'') subset of vertices 
Vo C V. Denote by Vj the vector of coordinates of vertex Vi G V. 

Definition 1.1 An equilibrium stress (or self-stress) is an assignment of real numbers Sij = Sji 
to the edges, a tension if the sign is positive, or a compression if the sign is negative, so that the 
equilibrium conditions 



hold at each vertex Vj e V\Vo (see Fig. 1.3). 




Figure 1.3 



A finite framework F{E; V) is called rigid in W'' if and only if there is a neighborhood N{F) 
of F in the space of parameters W''^ such that any other W'' realization of graph {E, V) from 
N{F) with the same lengths of all edges is congruent to F{E; V). Sometimes it is interesting to 
study rigidity of graphs some of whose vertices are pinned. Framework F(E; V, Vq) is called rigid 
with pinned vertices Vq if F has a neighborhood N{F) in M*^^ such that any other realization of 
{E; V) from N{F) with the same lengths of all edges and the set of pinned vertices inherited from 
the set of pinned vertices of F is congruent to F. If the "neighborhood condition" is dropped the 
framework is called globally rigid. A framework that is not rigid is called flexible. A framework 
F{E; V, Vq) that supports an all non-zero equilibrium tension is referred to as a spider web. 

Definition 1.2 An framework F{E; V, Vq) (possibly infinite) is referred to as rigid if any finite 
subgraph of F{E;V,Vo) is contained in some rigid finite subgraph of F{E;V,Vo). 

By polygonal partition of a planar set with piecewise-linear or no boundary we mean an 
edge-to-edge partition of this set into convex bounded polygons. The 1— skeleton of a partition 
is a framework whose vertex set is the vertex set of the partition, and whose edge set is the 
partition edge set. A polygonal partition which consists of triangles is called a triangulation. 




iij)GE 
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Definition 1.3 Let M be a set with polygonal boundary in (M might be all o/M^J, and let Ti 
be a collection of open polygons in M, such that the number of polygons intersecting any compact 
subset of W is fimte. We call the elements of % holes and denote by H the pointwise union of 
the holes. We say that M\H supports tension if M\H admits a partition with the set of edges 
E and vertices V , such that the framework {E, V,V D dM) is a spider web. 

Evidently, a polygonal partition can be replaced by a triangulation in this definition without 
loss of generality. A direct generalization of this definition to the case of general dimension is 
possible, but not quite natural, since not all spider webs in dimensions higher than 2 can be 
interpreted as 1-skeletons of polyhedral partitions (see Connelly and Whiteley (1993, 1996)). In 
the planar case the situation is simplified by the fact that any spider web with self-intersections 
can be turned into the 1-skeleton of a polygonal partition by adding points of self-intersections 
to the vertex set of the framework, and modifying the edge set accordingly. The cone of tension 
of the new partition contains the cone of tension of the original one. A somewhat more natural 
definition would the one in which the existence of some spider web in the complement of the 
holes is required. Let us now make some observations about holes. First, if a hole is non-convex, 
then there is no triangulation of the complement such that its 1-skeleton minus the boundary 
of M supports equilibrium tension. For instance, the equilibrium of forces at vertex v on Fig. 
1.4 is impossible if all edges incident to this vertex are under tension. Therefore, if two holes 
overlap, and their union is not convex, tension vanishes on all of their convex hull. This is called 
the coalescence effect of overlapping holes. 




Figure 1.4 

Notice that the above observation about the convexity of connected components of the com- 
plement of the set supporting tension is not valid for dimensions greater than three, whatever 
definition of tension we adopt. For instance, a three-dimensional polytonal hole can have a sad- 
dle point, to which a number of edges lying in the complement of the holes can be attached so 
that this point will be in static equilibrium. However, a set supporting tension in R'^cannot have 
points of strict convexity. 

It is natural to call subsets of M that cannot support tension defects. The primary defects 
are, obviously, holes themselves. The following definitions iteratively define the notion of a defect 
of /c*'* generation. Note that in these definitions the polygonality of the holes is not important. 

Definition 1.4 Let be a set of holes. Elements ofH are called defects o/O*'* generation. 
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Definition 1.5 A connectivity component (understood topologically) of defects of k*^ generation 
is referred to as a k— cluster. 

Definition 1.6 A defect of {k + 1)"' generation is the convex hull of a k— cluster. 




Figure 1.5 

Fig. 1.5 provides an illustration to our definitions. Here we used the solid line for the 
boundary of holes (= defects of O*'* generation), the dash line for the boundary of a defect of 
the first generation, and dotted line for the boundary of a defect of the second generation. Note 
that our definition of defect does not cover all subsets that cannot support tension, since an area 
that fails to support tension need not belong to the convex hulls of a connected component of 
the set of holes as it was noticed above (see Fig. 1.2 and explanations in Section 1.2). However, 
the effect of a "pinwheel configuration" is insignificant for the "naive" Poisson percolation model 
that we analyze in this paper. Even without taking this effect into account, we show that tension 
disappears on all of W . 

The existence of a triangulation which supports tension may have some interesting implica- 
tions for the modeling of physical properties of materials with networks of Hooke springs and 
geometry of convex surfaces (see Subsection 1.2). Among other things it implies that every tri- 
angulation of M\H is globally rigid not only in M^, but also in , provided all its vertices lying 
on the boundary of M are pinned (Connelly (1993)). 

The elasticity (rigidity) properties of a glass are related to how amenable the glass is to 
continuous deformations which require little energy. From a physical point of view, it is not 
enough to declare that the distance constraints force the structure to have only one configuration, 
since the bonds in a physical network do not behave as ideal edges in a framework. There should 
be a way of describing the behavior of the system as it is perturbed. That is why physicists 
often consider the energy function defined on the edges of a network of Hooke springs, i.e. each 
spring has some optimal length at which its energy is minimal, stretching or shortening a string 
increases the energy of this connection 

A tensegrity framework is a generalization of this model where besides Hooke springs there 
are members whose energy increases with the distance, and members whose energy decreases 
with the distance. In context of energy considerations it is often useful to work with the notion 
of tensegrity framework (Connelly, Whiteley (1996)). In a tensegrity framework all edges are 
partitioned into three types, cables struts and bars Eq, i.e. E = EqU E^ U E_. Together, 
struts, cables and bars are called members. Let F{E,Eq,Ej^,E_;V,Vq) be some tensegrity 
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framework in W'' , and denote by Ifj the length of a member (ij). If a cable is stretched the 

energy in the cable increases. If a strut is shortened, the energy in it decreases. Any change 

in the length of a bar gets the energy to increase. Therefore networks of Hooke springs are bar 

tensegrities from mathematical point of view. Thus the energy S^ij of member (ij) considered as 

the function of its squared length 

is monotone increasing if (ij) is a cable, 

is monotone decreasing if (ij) is a strut, 

has a strict local minimum at l^j called the equilibrium length. 

It is natural to define the energy function of a finite tensegrity framework as the sum of the 
energy functions of its members. Thus 



The simplest way to define the energy function when all members are bars is as follows 



where the sum is over all ordered pairs of vertices of the framework, riij is or 1 depending on 
whether there is a bond between i and j, Uj is the length of the bond between i and j, 1^^ is the 
equilibrium bond length, and a^j > is the spring constant of the bond (ij). 

In the spirit of the definition of equilibrium stress we assume that a strut can support only 
compression, a cable can support only tension, and a bar can be under both types of stress. For 
more detailed information on tensegrities see the works of Connelly and Whiteley (1995), and 
Connelly (1993). 

Consider now the energy ^ as a function of the coordinates of the vertices of the framework. 
A finite tensegrity framework F in M'^with pinned vertices Vq C V is called prestress stable if 

1) the differential of (considered as a function on the space of parameters K^V) at the 
point corresponding to F is zero, 

2) the second differential of ^ is a positive semidefinite quadratic form whose kernel re- 
stricted to infinitesimal motions leaving Vq unmoved consists of trivial infinitesimal motions of 
the framework. 

Definition 1.7 An infinite tensegrity framework is called prestress stable if every its finite sub- 
graph G is contained in a prestress stable subgraph whose un-pinned vertices contain the vertex 
set of G. 

The concept of prestress stability is due to Connelly and Whiteley (1993, 1996) and basically 
accounts for local minima of the energy function. It turns out that if the complement of H in M 
supports tensions, then for any triangulation T of M\H there is a an interpretation of the edges 
of T as struts and cable such that the resulting tensegrity is prestress stable in 3-space with dM 
pinned. In particular it implies that T is prestress stable as a bar framework in 3-space with dM 
pinned. 

So far we discussed only statics of framework. There is an interesting implication of the non- 
existence of tension in the complement of a collection of holes for mechanics. It was conjectured 




(1) 
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by Connelly (1998) that in this case the complement has a triangulation which is a flexible in 
3-space with boundary pinned. There also are interesting connections between our problem and 
convex geometry that in its original form are due to Maxwell (1864, 1869-1872) and Cremona 
(1872). They are outlined in Section 1.2. 

Let us summarize the implications of the existence of tension in the complement of the holes 
(see also Connelly, Mitchel, and Rybnikov (1998, 1999)). 

Proposition 1.8 Let M be a convex set of ^ with polygonal or no boundary, and let % be a 
collection of convex open polygons in M , possibly overlapping. Suppose M\{H U dM) supports 
tension. Let P be a convex polygonal partition of M\H with the set of vertices V and the set of 
edges E, and let F{E; V,V Cl dM) be the corresponding framework with pinned vertices V fl dM. 
Then 

1 ) framework F{E\ V, V fl dM) is globally rigid in as a bar framework; 

2) there is a partition of E into two sets E^ and E_, such that if members of E^ are cable 
and members of are struts, the resulting tensegrity is prestress stable in . 

Conjecture 1.9 (Connelly) Let M be a convex subset ofM? with polygonal or no boundary, and 
let Ti be a collection of convex open polygons in M. Suppose M\{H U dM) does not support 
tension. Then there is a triangulation of M\H which is flexible in 3-space with the boundary of 
M pinned. 

We say that tension vanishes at point p e M if there is no subset of M\H supporting tension 
that covers this point. In the following subsection we show that for a finite system V. any subset 
of M\H which supports tension is, in fact, contained in the maximum subset Smax ^ M\H 
supporting tension. Thus, when the number of holes is finite, M\H can be partitioned into two 
polygonal subsets, the maximum subset supporting tension and its complement where tension 
vanishes. We conjecture that this is also true for infinite case ( see Section 1.2). The above 
proposition suggests an interpretation of the area that supports tension as stable one, and of the 
area that cannot support tension as amenable to deformations requiring little energy. 

One of the applications we have in mind is percolation of defects in presence of crystallization. 
Here H. is not a set of holes but rather a dissemination of some alien substance that have poor 
contact with our material. Crystallization cannot take place in a region where the tension has 
been lost, but we assume that the system has some memory. The reason why we assume that 
crystallization cannot occur when tension has been lost is that the loss of tension means the 
loss of a rigid structure as it is suggested by Proposition 1.8 and Conjecture 1.9. We do not 
claim that the rigidity of flexibility of a triangulation of the complement of the defects directly 
explains the behavior of an actual physical system. However, if Connelly's Conjecture 1.9 is 
correct, the existence of tension is a necessary and sufp.cient condition of the global rigidity of all 
triangulations of a planar set with polygonal boundary. 

We assume that tension cannot be lost instantaneously, since there is some contact between 
the material and the alien substance. Before tension disappears on defects of /c*'* generation it 
disappears on defects of {k — I)*'* generation. In section 2 we show that there is N G Z+ such 
that with probability 1 every point of is covered by a defect of A^*'*^ generation. It can be 
interpreted as that there is time tc{7i) such that if the system has not crystallized by time tc(^) 
it will never crystallize. 
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1.2 Tension and convexity 

J. C. Maxwell (1864, 1869-1872) and later L. Cremona (1872) observed that a framework in 
which can be interpreted as the 1-skeleton of a 2-sphere admits an all non-zero equilibrium 
stress if and only if this framework is the projection of a polyhedral sphere in . This is also 
true for a partition of into convex polygons, but in this case an infinite polyhedral surface 
plays the role of a 2-sphere. These theorems were later rigorously proved by Crapo and Whiteley 
(1982, 1984, 1993). General theory of the relationships between stresses, liftings, and Voronoi 
diagrams can be found in Rybnikov (1999). In fact, this connection between stresses and convex 
surfaces provides much insight into the problem of understanding the redistribution of stress 
caused by punched holes. It turns out that in Maxwell construction a tension on the 1-skeleton 
of a polygonal partition corresponds to a convex lifting of the partition. Therefore all of the 
complement of the holes supports tension if and only if the holes can be lifted to different faces 
of some convex surface. Projecting such surface back on the plane we get a partition of the 
complement of the holes whose skeleton support an all non-zero tension. 

Proposition 1.10 Let H. be a discrete system of holes in W' . The complement of H supports a 
tension if and only if there is a continuous function defined whose graph is a convex polyhedral 
surface in such that all connected components of H are vertical projections of different facets 
of this surface 

For example, on 




Figure 1.6 
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For example, on Fig. 1.2 we see that no convex piecewise-linear function can lift the triangles 
to three different facets, since it would imply that these triangles can be mutually separated in 

by three lines meeting at one point. This is clearly impossible for the configuration depicted 
on Fig. 1.2. Similarly, on Fig. 1.6 we see that the complement of the original system of holes 
(thin solid line) does not support a tension. However, the system of enlarged holes (dash line) 
supports a tension. Holes A, B, and C form a "pinwheel configuration" like one on Fig. 1.2. 

As we have seen the Maxwell correspondence helps to understand the redistribution of tension 
in the case of a finite number of holes. In fact it also suggest an algorithm for computing the 
maximum subset of the complement that supports a tension (Connelly, Mitchel, Rybnikov (1998, 
1999)). The following theorem shows that one can exactly define the maximum subset of the 
complement that supports a tension, and illustrates the use of Maxwell correspondence. 

Lemma 1.11 Let M be a convex polygon in (M might be all of^), and let V. be a finite 
collection of convex polygons in M. There is a set Hmin such that 

1) H C Hmin, 

2) M\Hjnin supports tension, 

3) There is no proper subset of Hmin satisfying the first two conditions. 
Furthermore, Hmin is unique. 

Proof. Let H' be a set of holes that satisfies the first two conditions, namely, H C H' and M\H' 
supports tension. Recall, that any system of holes satisfying the first two conditions is actually 
a collection of edge-disjoint convex polygons, i.e. no pair of polygons in the collection share an 
edge. 

Suppose not all of the vertices of H' come from the vertices of H. We show that in this case 
H' is not minimal, which means that there is a subset of H' that satisfies Conditions 1 and 2 and 
is contained in H'. Let -D be a polygonal hole of H' with a vertex v which does not belong to the 
vertex set of H. Vertex v can be separated from all other vertices of D by a line passing though 
some two vertices of H. The resulting triangle A can be cut off D without violating the property 
of M\{H'\A.) to support tension. Let us illustrate it with the use of Maxwell correspondence. 
M\H' has a triangulation that lifts to a convex surface. Denote the lifting map by L. Let us 
modify the surface by cutting off the corner of the vertex L{v) by a plane passing through the 
other two vertices of A and cutting through all other edges incident to L{v) at a sufficiently 
small distance from L{v) (see Fig. 1.7). 
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Figure 1.7 

By Maxwell correspondence the projection of the resulting surface back on the plane gives 
a triangulation that supports tension in its 1-skeleton. Thus, if minimum subset exists then it 
consists of holes all whose vertices come from the vertices of H. 

There is a finite number of way to construct a system of holes in M that cover H, whose 
vertices come from vertices of H and M, and whose complement supports tension. Pick among 
such systems of holes a system with minimal area. Denote it by Hi. We claim that this is, indeed, 
Hmin- Assume that there is a proper subset H2 of Hi that satisfies the first two conditions. 

The procedure described above can be repeatedly applied to H2 until the vertex set of what 
is left of H2 is a subset of the vertex set of H. But then the resulting collection of holes satisfies 
the first two conditions, and has smaller area than Hi. It contradict to the choice of Hi. Thus 
Hi is minimal with respect to inclusion. Now, let's show that Hi is actually contained in any 
system of holes satisfying the first two conditions. Suppose H' is another system of holes which 
does not contain Hi, but satisfies the first two conditions. Consider a system of holes whose 
holes are intersections of the holes of H' and Hi. Denote it by Hmesh- Since both H' and Hi 
covered H, Hmesh covers H too. 

By definition M\Hi and M\H' both support tension, i.e. have partitions whose 1-skeletons 
support tension. Denote these graphs Gi and G'. Consider a new graph, whose vertex set consists 
of all the vertices of the old graphs, Gi and G', and all intersections of edges of the old graphs 
( this graph is what we see when Gi and G' are drawn on the same picture of M). The edge 
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set of the new graph is defined accordingly. Denote it by Gmesh- Tensions on Gi and G' give 
rise to a tension on Gmesh- Thus, Gmesh gives a partition of M\Hmesh whose 1-skeleton supports 
an all non-zero tension. Therefore we have a system of holes on M that satisfies the first two 
conditions, and is contained in Hi. But we have already shown that this is impossible. Thus 

Hi = Hmin- ^ 

If H is infinite we can still conclude that a maximal (with respect to inclusion) set supporting 
tension is the complement of a set of edge-disjoint polygons all whose vertices come from the 
vertex set of H, since the argument about cutting a corner can be also applied to an infinite 
surface. 

As a consequence of the above theorem, the problem of finding the maximal area supporting 
tension can be rephrased as the following optimization problem on convex polyhedra. Find a 
convex surface in such that each hole of T-L is covered by the vertical projection of a facet of 
the visible from the plane part of this surface, and the total area of the projections of the facets 
that cover holes is minimal possible. Another important implication of this theorem is that for 
finite systems of holes the complement of the set where tension vanishes must support tension, 
i.e. it has a triangulation which is a spider web. The case of infinite system of holes is more 
complicated. Even under additional restrictions on the system of holes, like that the vertices of 
the holes form an (r, i?)-system, or that the sizes of the holes are uniformly bounded from both 
above and below it is not obvious that the union of all subsets of holes supporting tension can 
be represented as the complement of a discrete set of edge-disjoint polygons. 

Conjecture 1.12 Let % he an infinite discrete system of polygons in W' . Then the union of 
all subsets of^\H supporting tension can be represented as the complement of a discrete set of 
edge-disjoint convex polygons. 

In addition we conjecture that the complement of the set where tension vanishes is not only 
polygonal, but also supports tension. Notice that our main result is about the vanishing of 
tension and the truth or false of this conjecture are irrelevant to our proof. 

2 Tension Percolation Models 

When not all of the complement of H supports tension, it is natural to ask what is the maximum 
subset of S\H supporting tension. 

Connelly and Rybnikov (1998, 1999) showed that the problem of finding the subset of max- 
imum area supporting tension can be reduced to a linear programming problem, and as a con- 
sequence there is a polynomial (in the number of vertices of the holes) algorithm finding the 
maximal tensed subset. This polynomial algorithm can be effectively implemented and used for 
computer simulation of the membrane model for different distributions of holes and their param- 
eters. However, this empirical evaluation for finite systems does not allow us to understand the 
limiting behavior of the system when the number of holes is infinite. In this section we consider 
a number of probability models for the distribution of holes and their sizes. 

In the previous sections we looked at the behavior of tension on finite membrane with polyg- 
onal defects. A natural question is how the tension is redistributed on an infinite membrane 
(film) with polygonal holes. One can assume different distributions of the holes and their shapes. 
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The naive model would be that all defects are centered at nodes distributed according to Poisson 
law, and the size of the hole is independent of the Poisson distribution. For this model we give 
a complete analysis and show that assuming that the holes are distributed according to Pois- 
son law with parameter A, and the shapes of the holes are IID and independent of the Poisson 
distribution, the tension ceases to exists on all of for any value of Poisson parameter A. 

A more realistic approach would be as follows. Consider a Voronoi partition of the plane 
whose set of sites is a realization of a Poisson distribution with parameter A. Declare each 
Voronoi domain normal with probability p and defective with probability 1 — p (for percolation 
on Voronoi partitions see Vahidi-Asl and Wierman (1990)). For each defective domain define 
a convex hole centered at the Voronoi center of this domain. We suggest a family of model in 
which each convex polygonal hole is situated at a small neighborhood of its Voronoi domain. We 
associate with each hole three distributions: the distribution of the number of vertices of the 
hole, the distribution of the angular coordinate of a vertex of the hole, and that of the radial 
coordinate of the vertex. Let z/ be any discrete distribution for the number of vertices of a defect. 
Let p be a distribution for the length of the radius-vector of vertex. A reasonable choice for p 
would be some random variable taking on all values between the radius of the inscribed circle 
and the radius of circumscribing circle of the domain. A natural choice for the distribution of 
the angular coordinate is a uniform distribution on the circle. The basic question is whether 
there is a critical value of parameter p such that there is an infinite connected component of 
supporting tension. A special case of such model would be one where the holes coincide with 
their Voronoi domains. We conjecture that in this case for any value of A and for any value 
of p there is no infinite connected component supporting tension. The constructions involving 
Voronoi partitions are motivated by our desire to have a model in which the effect of pinwheel 
configuration is significant. In the "naive" Poisson model the vanishing of tension is caused by 
the overlapping of holes, which is in turn is caused by the independence of the distribution of 
the hole sizes and the underlining Poisson process. 

A lattice analog of continuous tension percolation models was introduced by Connelly and 
Rybnikov (1998, 1999). Let T be a triangular lattice on the plane. Remove each edge with 
probability 1 — p. The question is whether there is a critical value of probability such when if 
p is greater than there is an infinite subgraph of T supporting tension. This looks reasonable for 
modeling tension in a membrane, but the assumption that all atoms lie on three families of parallel 
lines does not seem to be realistic for describing properties of molecular systems like glasses. Thus, 
a generic version of this problem should have instead of T a generic triangular lattice. Note that 
the problem of tension percolation on a lattice is related to rigidity percolation on lattices. 
See Holroyd (1998) for rigorous estimates of the critical probability of rigidity percolation and 
Duxbury, Jacobs, Thorpe, Moukarzel (1995, 1996, 1997, 1998, 1999) for simulation results and 
their physical interpretation. 

3 Percolation on iterative defects 

In this section we show that if the centers of holes are distributed in according to Poisson 
law and their shapes are IID, tension disappears on all of M^. Moreover, the following results 
on iterative convex hulls are true not only for dimension two, but also for W^. The proofs can 
be adopted straightforwardly. It is convenient for as to adopt a more general definition of hole. 
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namely without the polygonality restriction. 

Definition 3.1 A hole (f-hole) centered atp^W'' is a region 

^^(p,/) = {p + /(p)x I ||x||<l} 

where f is a continuous positive function defined on a unit {d — l)-sphere. 

Let V. be some collection of holes. 

Definition 3.2 We say that a collection of holes Ti = {Hi{pi, /i), i?2(p2, f2), ■ ■ ■} is uniformly 
bounded from below, if there exist r > such that /j(x) > r for all fi and ||x|| = 1. Similarly, we 
say that % is bounded from above, if there exist R> such that /i(x) < r for all fi and ||x|| = 1. 

Consider a d— dimensional Poisson point process with rate A. Let Y = Y(uj) be the collection 
of nodes of some realization co of the process. Each node y e y{co) is the center of a hole H{y, fy), 
where function fy is positive and continuous. 

Let n he a probability measure on some subspace of C(S''^^, R_|_). Suppose that for each 
y the function fy is chosen from a distribution // independently of the other functions and the 
configuration u. Therefore, the holes H(y, fy) are IID. Under this condition, it is known that if 
v. is bounded both from above and below, there exists some Ac = Ac(?^, //) such that for A < Ac 
there exists no infinite cluster and for A > Ac such cluster exists a.s (see Menshikov and Sidorenko 
(1987)). 

The next statement is our main percolation result. In part, it can be viewed as a generalization 
of some results of Schonmann (1992) for bootstrap percolation on Z*^. 

Theorem 3.3 For any distribution /j, and any A > there exists a nonnegative integer N = 
N{ii, A) such that a.s. there is an infinite N -cluster whose convex hull coincides with . 

To prove this theorem, we reduce it in two steps (see Lemmas 3.4 and 3.5) to a simpler 
proposition (see Proposition 3.7). 

Lemma 3.4 It is sufficient to prove Theorem 3.3 only for the case when the set of defects is 
bounded from below by some r > 0. 

Proof. According to definition, a hole corresponds to a continuous 

and positive function f{u) on a unit sphere S'^^^. Since S'^^^ is compact and f{u) is con- 
tinuous, f{u) achieves its minimum there, and this minimum is strictly positive. Let p{f) := 
min„g[//('u) > 0. The distribution on functions / induces a distribution on a real- valued 
positive random variable p. Therefore, there must exist r > such that P{p > r) = po > 0. 

Now we couple the original process with rate A and the functions fy distributed according 
to p, with a process which culls all the defects for which p(/) < r. It is straightforward that 
this process is also a Poisson process, but with rate A = Xpo. The distribution of defects p for 
the thinned process is such that the radius of the inscribed circle for any defect is at least r. 
Consequently, the set of holes is now uniformly bounded below. If we show that the statement 
of Theorem 3.3 holds for p and any A > 0, it will imply Theorem 3.3 for a generic p, as the set 
of holes in "tilde" model is stochastically smaller than that of the original one. □ 
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Hence, we can assume that each hole (= defect of 0*'* generation) corresponding to a node 
y & Y contains a ball of radius r with the center at y. Since taking a convex hull is a monotonous 
operation (namely, if A C B then conv{A) C conv{B)), it suffices to prove Theorem 3.3 for the 
case when all holes are the balls of radius r. Moreover, without loss of generality we will can 
assume that r = 1 as we can always re-scale the space. Consequently, Theorem 3.3 follows from 
the next lemma. 

Lemma 3.5 Let H, be a collection of balls of unit radius whose centers are distributed in 'M^ 
according to a Poisson law with parameter A > 0. Then there exists a nonnegative integer 
N = N{\) such that a.s. there exists an infinite N -cluster whose convex hull is . 

3.1 Proof of Lemma 3.5 

We will present the proof for the case c? = 2 (so that the defects are in K-^). The arguments can 
be generalized for d > 3 rather easily. 

Pick a hole(which is a unit circle) from 1-L, and denote it by G(0). Examine all circles from 
% which have non-empty intersection with G(0), take the convex hull of them and (5(0), and 
denote the resulting set by G{1). Then consider all circles of l-L which intersect G(l), and call 
the convex hull of G{1) and their imion G(2). Iterating this procedure we construct the sets 
G(2), G(3), .... If for some k no element of % intersects with G{k), we have G{k + 1) = G{k) for 
all I > 0. 

Proposition 3.6 G{k) is contained in some defect of k^'^ generation. 

This statement is obvious, because the operation of taking convex hull is commutative. 

For any point a e let {pa,'^a) be its polar and {xa,ya) be its Cartesian coordinates. 
Without loss of generality assume that G{0) is centered at the origin. Denote by C{k) the circle 
of radius k centered at the origin. 

Proposition 3.7 With a positive probability depending on X only, 



that is, G{k) 's eventually cover all as k ^ cxd. 

Proof. Here we actually prove a stronger result: we show that (3) holds even if G{k + 1) were 
constructed using G{k) and only those circles of T-L whose centers lie inside the ring 



G{-) is an increasing process taking values in the subsets of ^ . Observe that G{k — 1) C 
C{k + 1)), and let be the event X{k + 1) C G{ky . Observe that the "smallest" value G{k) 
can take on is C{k + 1). Therefore, if we can demonstrate that 



C{k + 1) C G{k) for all A; = 1, 2, 3, . . . 



(3) 



Rk+l 




P{Ek+i I G{k) = C{k + l),G{k-l),G{k-2),..., G{1)) 
= P{Ek+i\G{k)=C{k + l))>l-jk, 



(4) 



15 



it will imply that P{Ek+i \ E^, E^^i, . . . ,Ei) > 1 — 'jk and as a result 

oo oo 

P{f]E,) = llil-^,)>0 (5) 

as soon as 7;^ < oc. 

It is straightforward that P(i?fc+i | ^(A,') = C{k + 1)) > for all k, so it suffices to estimate 
from (4) only for large values of k. Suppose that G{k) = C{k + 1) and break the ring into 
M = M{k) := 2'K\lk + 3 congruent pieces 

4+1 = {ciei?;t+i : e (m,a(i + l)]}, i = 0, 1, . . . , M - 1 
where a. = a{k) 



M V^T3 

Let the event E^^-^ be contains at least one center of a circle from . Then these events 

are independent for different values of i and k, and consequently the probability of the event 
Ek^i = ^-^i+i "each -R^*|i contains such center" is 

P(^,+i) = (1 - e-^"('=+7/4)/2)^ > 1 _ Me-""('=+^/4)/2 ^ ^ _ c^j,il2^-c,k^'- 

for some positive constants Ci and C2 and large k. 

Our next step is to show that E^+i implies E^^i whenever C{k + 1) C G{k). Once this is 
established, we have 



Ik 



< c,k'i\-'^^'"\ 



and, as a result, ^7^; < 00, which yields (0) and proves the proposition. 
Assume that E^+i occurs and consider an arbitrary point 

aeC(A; + 2). (6) 

Since each of contains a center of a circle of H and every such circle obviously intersects 
C{k + 1), there exist two points b and c in G{k + 1) such that 

Pb>k + 5/2, pc>k + 5/2, 
^Pb<^a< mod 27r , 
'Pa ~ fb ^ oi mod 27r , (Pc ~ 'Pa ^ oi mod 27r . (7) 

Since G{k + 1) is convex, it must contain all points inside the triangle formed by the origin 
and points b and c. For any two points e and / of the plane, let [e, /] denote the segment with 
these points being its endpoints. There exist a point h lying on the segment [b, c] such that [0, h] 
is orthogonal to [b, c] . Notice that this point has the smallest distance from the origin among all 
points of the segment [b, c]. Let us get the lower bound on p/j. From (7) it follows that at least 
one of the angles between [0, h] and [0, b] or [0, c] does not exceed a. Suppose it is the first one. 
The length of [h, b] in this case is smaller or equal to 

^ Pb 
pb sm a < phOL 
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whence 



as (A; + 5/2)^ > {k + 2){k + 3). Consequently, equations (7) and 6 yield that the point a lies 
inside the triangle A(Wc) and hence belongs to G{k + 1). □ 

Let us make some observations which rely on obvious generalizations of Proposition 3.7. In 
our discussion we will rely on the proof presented above. 

First, color the plane as an infinite chess-board with the cell size of 1/4, so that a point with 
coordinates {x, y) is green if [4xJ + [4|/J is even, and blue if this sum is odd. The arguments in 
the proof will remain valid if instead of Ti we consider only those circles of it which have "green" 
centers (these centers form a Poisson point process with the same rate A on the "green part" of 
the plane). The similar statement about the "blue" process is valid as well. 

Next, note that if while constructing G{k + 1) from G{k) we do not use any circle with a 
center a such that 

-A;-^/^ <<^a< k'"^'^ mod 2ti or 

TT - A;-2/3 < ^„ < TT + A;-2/3 (8) 

(and use only the green part of the plane), the arguments above still hold for large A;, as A;~^/^ = 
o{a{k)). Denote the resulting process G'{k), and the probabilities corresponding to (4) as 1 — 7^. 
Since 7^ are asymptotically of the same form as we have ^7^ < 00. 

The third observation is that if the circle C{ko) {ko > 1) is completely covered by some 
elements of ^ with all their centers lying inside C{ko), then the probability p(A:o) that the process 
G'{k) described above (using only green part of the plane and "avoiding" points close to the 
horizontal axis), can be continued indefinitely, in a way described by (3), has the property 

lim p{ko) = 1. 

kg— ^00 

This follows from the fact that p{ko) = Yi'kLkoi^ ~ Tifc) ^ 1 ~ '^k>ko ^'k right 
hand side converges to zero as A;o ^ 00. 

Consequently, there exist ko > 1 such that p{ko) > 0.9. Fix ko and let v denote the probability 
that indeed C{ko) is completely covered by circles of?/ as described in the previous paragraph. 
The value 1/ may be small, however it is positive and therefore there exist a positive integer T 
such that 

(1 - uf+' < 0.1 (9) 

Consider T + 1 circles of radius ko on the plane with their centers located at the points 

(0,0), (A;o^O), i2kl0), (TklO). (10) 

According to (9), the probability that at least one of them is completely covered by circles of 
v. is larger than 0.9. Arbitrarily pick one such circle (let it be, say, the i"' one). The event 
"the process G'{k) starting from {ikQ,0) can be continued indefinitely" is independent of the 
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configuration inside the other circles, since when the process G'(-) reaches the j*'* circle {j / i) 
k ^ \i — j\kl =: mkl, m > 1, and therefore the j*'* circle lies entirely inside the angle 

111 1 

m/cQ ^ mkQ mk'^ ^ mkQ 

and at the same time 

1 1 1 

< ~ 

mkl ~ m?l^kl k'^l^' 

Hence, with probability exceeding (1 — (1 — i^)^"*"^) x p(A;o) > 0.81 one of the A;o-radius circles 
with a center at one of the points of (10) is completely covered by elements of ?^ anrfthe process 
G'{k) which starts from z"' circle can grow indefinitely. 

Now let us tile the plane with the boxes L x L where L = STfeg, so that for any G Z^, 
the box B{i,j) consists of the points 

{ {iL + x'jL + y') where {x' , y') e (0, Lf } 

Suppose that i + j is an even number. We say that B{i,j) is open (see Fig. 3.8), if 

(A) one of r + 1 circles of radius ko with the centers in the set 

{Xi, Yj), (Xi + kl Yj), (Xi + 2kl Yj), (X, + Tkl Y,) 

^{{t + 0.6)L,Yj) 

where Xi = {i + OA)L, Yj = {j + 0.5)L 
is completely covered by the elements of V. and 

(B) the process G'{-) starting with the center aij of a circle satisfying (A) and using only 
green points for k > ko reaches the circumference of radius of 0.6L (i.e. G'{ki — 1), where 
ki = [0.6L], contains the entire circle a^j +C{ki) ). 



green 


blue ^.-'^ 


green 


blue 


green / 




\ 
/ 


blue 


green 


L 




blue 


green 



R=0.6L 
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Figure 3.8 



In the case when i + j is odd, we call B{i,j) open if the same is true except G'{k) uses blue 
points. 

From the translation invariance and the previous arguments it follows that 

P{B{i,j) is open) > 0.8 

Moreover, the boxes B{i,j) are open independently of each other, since the neighboring (by an 
edge) boxes use points of different colors and G'{ki — 1) does not reach any of the corners of the 
box (since 0.6L < y^(0.5L)2 _|_ (o.4L)2). Also observe that when two boxes B{i,j) and B{i',j') 
sharing a common edge are open, the circles C{ki) + a^j and C{ki) + ai>j> have a non-empty 
intersection. 

We say that two boxes are connected, if there exists a sequence of open boxes starting with 
the first one and ending with the second one such that each pair of two consecutive boxes in 
this sequence share an edge. The "box"-cluster containing the box B{i,j) is the set of all boxes 
connected to it. 

Hence, we can couple the boxes {B{i,j)} with the site percolation process on where each 
site is open independently of the others with probability exceeding 0.8. However, as follows from 
Zuev (1987), the critical probability Pc for the site percolation on is smaller than 0.68 < 0.8 
and therefore there exists a.s. an infinite open cluster of boxes B{i, j). As a result. Proposition 3.6 
yields that there exist a.s. an infinite {ki — l)-cluster. 

Let us prove the following 

Proposition 3.8 Consider a site percolation model on 1? in the supercritical regime (i.e. p > 
Pc). Then the convex hull of a (unique) infinite cluster is M^. 

Proof of the proposition. First, we remark that the fact that the infinite cluster is unique whenever 
it exists follows from Aizenman, Kesten and Newman (1987). Let C(oo) be such infinite cluster. 
On the plane with coordinates {x, y) consider the angle 

A = {{x,y) eZ"^ : x>0, 0<y <x}. 

Notice that the plane can be partitioned into 8 cones congruent to A. Let Ea be the event 
|AP|C(oo)| = oo, i.e. there are infinitely many points of the infinite cluster lying inside of A. 
Since this is a tail event, P{Ea) is either or 1. If it were 0, it would imply that the probability 
that there are infinitely many points of C(oo) in each of the other seven congruent to A cones is 
also 0. This contradicts to the fact that |C(oo)| = oo. Hence, P{Ea) = 1 and the similar is true 
for the others cones. However, it is easy to see from geometrical observations that if each of this 
eight cones contains infinitely many points of C(oo) then conv{C{oo)) = E? . □ 

Hence, the convex hull of the {ki — l)-cluster is at least as large as that of an infinite open 
cluster for the site percolation on Z^ (rescaled L times), and Lemma 3.5 has been proven. □ 
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